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BLAST  PHENOMENA  FROM  EXPLOSIONS 
AT  THE  WATER  SURFACE 


A.  SAKURAI,  J.  M.  PINKSTON,  JR.,*  AND  J.  N.  STRANGE 
NUCLEAR  WEAPONS  EFFECTS  DIVISION  1 

U.  S.  ARMY  ENGINEER  WATERWAYS  EXPERIMENT  STATION 
VICKSBURG,  MISSISSIPPI 


Introduction 

'  The  principal  objective  of  this  study  was  to  evaluate  the 
blast  and  shock  phenomena  (in  air  and  underwater)  resulting  from  a 
pure  surface  explosion  (i.e.  the  situation  where  the  charge  is  half 
in  and  half  out  of  the  water). 

The  particular  case  of  an  explosion  at  an  air-liquid  (water, 
in  most  cases)  interface  is  considered  in  this  paper  because  this 
case  is  easier  to  handle  both  analytically  and  experimentally  than 
an  explosion  at  an  air-solid  interface,  and  thus  permits  one  to  com¬ 
pare  more  readily-  analytical  and  experimental  results.  Moreover, 
it  should  be  noted  that  solids  under  sufficiently  strong  impulsive 
loads  are  known  to  behave  much  as  liquids.  The  general  results  of 
the  air-liquid  interface  case  can  be  utilized  for  any  air-solid  case 
by  substituting  the  appropriate  material  constants  to  represent  the 
liquid  properties  of  the  solid  under  consideration.  The  problem  of 
predicting  surface  burst  phenomena  involves  particular  difficulties 
when  compared  with  the  so-called  free-space  cases  in  which  the  explo¬ 
sion  occurs  entirely  within  a  given  surrounding  medium.  The  most 
serious  difficulty  arises  from  the  fact  that  the  general  character¬ 
istics  of  the  shock  waves  vary  depending  or.  the  nature  or  the  type 
of  the  explosion,  since  the  majority  of  the  features  of  the  shock 
waves  are  controlled  by  the  initial  stages ,  where  the  characteristics 
of  the  individual  explosion  predominate.  This  situation  makes  it 
difficult  to  establish  any  simple  scaling  rule  independent  of  the 
individual  nature  of  the  explosives  used.  For  this  reason  it  is 
usually  difficult  to  deduce  anything  ol  a  general  nature  from  the 
results  of  individual  tests  or  step-by-step  numerical  solutions 
starting  from  specified  initial  input  conditions.  A  less  serious 
difficulty  is  found  in  the  magnitude  of  the  density  ratio  of  the  two 
interface  media  involved.  Although  this  ratio  for  air  to  water  is 
much  less  than  1,  it  is  not  small  enough  to  warrant  neglecting  it. 
Thus,  the  motions  in  the  different  media  cannot  be  treated  separate¬ 
ly,  as  in  the  case  of  air  to  ground. 


vtw  u  »+•.  y:  v~'  -  * £ 5 ^ 


SAKURAI,  PINKSTON,  AND  STRANGE 

Aside  from  existing  literature  on  the  air-ground  cases ,  only 
a  few  previous  studies  on  the  air-water  case  are  available  for  either 
experimental  or  theoretical  investigations  (references  1  through  9) • 
All  these  data  are  based  on  rather  small  chemical  explosions ,  and  the 
relation  of  these  results  to  nuclear  burst  cases  is  very  uncertain, 
especially  for  this  particular  configuration  (i.e.  a  spherical  charge 
half  submerged  in  water) . 

This  situation  thus  implies  a  definite  need  for  a  thorough 
study  of  the  problem,  especially  for  nuclear  burst  cases.  It  is  es¬ 
sential  to  study  the  initial  and  later  stages  separately,  so  that  the 
proper  approaches  can  be  sought  for  handling  each  of  these  different 
stages,  whose  characteristic  features  are  depicted  in  fig.  1. 

As  mentioned  earlier,  there  may  be  no  general  solution  that 
will  represent  all  the  initial  stages  caused  by  different  kinds  of 
explosions,  such  as  chemical,  nuclear,  electric  spark,  etc.,  so  that 
they  must  be  studied  individually.  On  the  other  hand,  some  of  the 
important  characteristics  of  the  flow  fields  in  the  later  stages  can 
be  expressed  in  a  general  term  regardless  of  the  nature  of  the  ex¬ 
plosions  ,  mainly  because  only  a  small  interaction  exists  between  the 
flow  fields  in  water  and  those  in  air  regions  above  the  water;  thus 
they  can  be  formulated  separately  in  the  following  manner.  A  linear 
acoustic  approximation  is  used  to  satisfy  the  water  environment, 
while  a  modified  free  airblast  approximation  is  used  to  express  the 
flow  field  in  air.  This  model  leads  to  a  set  of  air  and  water  shock 
pressure  formulas ,  which  contain  some  unknown  parameters  and  func¬ 
tions  to  be  determined  from  the  characteristics  of  the  initial 
stages,  depending  on  the  type  of  explosive. 

Two  di/ferent  types  of  initial  stages  were  considered: 
a  point  source  model  of  a  finite  amount  of  mechanical  energy  released 
instantaneously  at  a  point  on  the  interface  assuming  this  to  be  real¬ 
istic  for  the  vide  varieties  of  nuclear  explosions ;  and  a  general 
chemical  explosion  of  bare,  spherical  HE, 

These  two  cases  were  utilized  for  determination  of  the  un¬ 
known  parameters  and  functions  appearing  in  the  pressure  formulas  in 
the  later  stages.  The  pressure  values  given  by  the  formulas  thus 
obtained  for  the  HE  cases  were  compared  with  the  test  data  from  the 
explosion  of  forty  10-pound  spherical  TNT  charges  conducted  simulta¬ 
neously  with  the  theoretical  investigation.  The  comparisons  of  the 
theoretical  values  with  experimental  data  and  with  test  data  avail¬ 
able  from  charge  weights  as  high  as  10 ,000  pounds  show  good 
agreement . 


Solution  of  Later  Stage 


General  characteristics  of 

the  flow  based  on  acoustic  theory 

Assume  here  linear  acoustic  properties  throughout  the  entire 
space.  This  assumption  is  obviously  not  valid  near  the  explosion 
source,  especially  in  the  air  environment.  Nevertheless  this  serves 
to  give  an  oVerall  insight  into  the  understandings  of  the  problem 


3 


SAKURAI,  PINKSTON,  AND  STRANGE 


which  can  thus  lead  to  the  development  of  modifications  necessary  to 
obtain  more  realistic  solutions. 

As  illustrated  in  fig.  2,  designate  p(r,z,t)  and  p'(r,z,t) 
as  the  water  and  air  overpressures  at  the  time  t  .  Let  r  and  z 
be  the  cylindrical  coordinates;  pQ  ,  p0  ,  and  cQ  are  the  pressure, 
the  density,  and  the  sound  velocity  in  the  undisturbed  water  region, 
and  pQ  ,  £T0  ,  and  TFq  Eire  the  corresponding  values  in  air. 

General  solutions  for  p  and  p’  ,  expressed  in  the  Laplace- 
Bessel  integral  forms,  are  fitted  together  at  the  interface  to  sat¬ 
isfy  the  boundary  conditions,  and  then  are  simplified,  by  utilizing 
the  small  value  of  the  density  ratio  e2  (e2  =  "tPo/Po)  >  t0  yield 

p(r,z,t)  =  F0(t  ~  ^H^t  -  -y  +  pa  (1) 

p’(r,z,t)  =  J  P0WQ[t  -  §-j H^t  -  y  (2) 

where  q  =  *Vr^  +  z^  ,  FQ  and  WQ  are  the  input  functions  related  to 
the  singularities  of  p  -  p'  and  3/3t(v'  -  v)  at  the  center.  Here 
v  and  v'  are  the  z  components  of  the  velocities  in  the  water  and 
air  regions,  respectively.  The  input  functions  F0  and  W0  are 
given  as 

P0(t)  sjij  (p’  -  p)2=0  rdr 

Wo(t)  It  <T’  -  ->K=0  rdr 

H  is  the  Heaviside  step  function,  and  pa  represents  the  term  cor¬ 
responding  to  the  water  pressure  induced  by  the  air  pressure  p* 
given  by  Equation  2.  Equations  1  and  2  exhibit  the  following  impor¬ 
tant  general  characteristics  of  the  flow  field.  First,  the  airblast 
field  p'  is  represented  by  a  spherical  wave  from  a  point  source; 
second,  the  water  shock  p  consists  of  two  parts:  one  caused  by  a 
doublet  source  and  the  second  one,  pa  ,  induced  by  the  air  pressure. 
While  these  qualitative  descriptions  of  the  flow  field  are  generally 
true,  since  they  are  based  primarily  on  the  small  value  of  e  ,  it 
is  only  necessary  to  modify  Equations  1  and  2  so  that  the  resulting 
formulas  can  be  used. 

Modification  of  the  formulas 

First,  Equation  2  for  the  air  pressure  p'  should  be  re¬ 
placed  by  the  equivalent  free  airblast  characteristics  caused  by  an 
explosion  in  free  air  with  an  equivalent  charge  weight  Wgf  which 
is  larger  than  the  actual  charge  weight  W  ,  and  the  ratio  6 
defined  by 
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Here  8  is  a  parameter  that  depends  on  the  nature  of  the  explosion. 
Second,  the  airblast  field  is  not  exactly  in  spherical  symmetry  but 
should  be  distorted  due  to  the  interaction  effect  with  the  water 
shock  field.  Since  this  distortion  is  rather  small,  this  can  be  in¬ 
corporated  in  the  formula  by  introducing  a  parameter  S  ,  which  will 
be  called  the  "Skewness  factor."  With  use  of  these  parameters  3 
and  S  ,  the  characteristics  of  the  air  pressure  p'  such  as  its 
peak  pressure  P(r,z)  and  its  duration  D(r,z)  at  (r,z)  are  deter¬ 
mined  from  the  corresponding  free-air  peak  pressure  Pf(Ae)  and  the 
duration  D^Ag)  at  an  equivalent  reduced  distance  Ae  as 


D(r,z)  =  8Df(Ae) 


(5) 

(6) 


where 


tan  0  *  -  (o  <  0  < 

z  \  -  -2} 


=  BqW 


1 

3  HS 
1  +  S  cos  0 


Using  this  -airblast  pressure  on  the  interface  as  an  input, 
the  air-induced  pressure  term  pa  can  be  estimated  numerically  by 
utilizing  the  scheme  given  in  reference  10;  here  alternatively  the 
following  approximation  formula  is  developed  by  introducing  various 
minor  assumptions 


t,a<r>z-t)  = 

o  o 


f  (f)  W*  - {>K<‘  - 


where  UQ  is  the  initial  airblast  velocity,  m^  is  its  initial  decay 
rate,  v  is  a  parameter  whose  value  is  about  0.1,  and  t(=t(r,z)) 
is  the  arrival  time  of  the  wave  at  a  position  (r,z).  The  resulting 
pressure  p(r,z,t)  consists  of  two  waves  represented  by  the  two 
terms  in  Equation  1  complemented  by  Equation  8.  In  principle,  the 
second  wave  induced  by  the  airblast  propagates  ahead  of  the  direct 
wave  of  the  first  term.  But  the  difference  is  very  small  for  most 
of  the  cases  (is  zero  for  r  =  0)  and  can  be  noticeable  only  in  the 
region  close  to  the  surface,  where  the  pressure-time  histories  of 


test  data  indeed  show  two  peaks.  It  is  also  noticed  that  the  first 
term  is  proportional  to  z  •  q~3  while  the  second  is  roughly  nro- 
portional  to  q-1  .  Thus  the  pressure  is  dominated  by  the  second 
term  with  increasing  distance  q  .  But  the  retaining  of  the  first 
term  is  indispensable  to  represent  the  importance  of  the  direct 
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effects  near  the  source.  Now,  the  input  functions  F0(t)  and  WQ(t) 
can  he  estimated  approximately  as 


F  (t)  =  (P  R  -  P’R'l 

0  l  s  s  s  s/r= o 

p"  W  (t)  =  (m'p'  -  e^mP  ) 

o  °  \  s  s/  , 


where  Ps  and  R_  are  the  water  shock  front  pressure  and  the  shock 
front  position,  P^  and  are  the  corresponding  values  for  the 
air  shock  front,  m  and  m'  are  the  decay  rates  of  the  water  and 
air  shock  velocities ,  and  a  is  the  effective  radius  of  the  spread 
of  these  inputs. 

Further  details  of  the  pressure  formulas  depend  on  the 
nature  of  an  explosion  through  the  determination  of  those  constants 
and  functions  introduced  above. 


Initial  Stage  and  Pressure  Formulas 


Bare,  spherical 
chemical  explosive 


Consider  here  the  case  of  the  spherical  chemical  explosive 
whose  center  is  positioned  at  the  water  surface.  The  magnitude  of 
a  in  this  case  can  simply  he 


a  =  aQ  =  "the  charge  radius" 


and  one  may  postulate  that 


Now  those  pressure  values  at  r  =  aQ  can  he  computed  using  the  free 
airhlast  and  free-water  shock  data  in  the  following  manner 


>  (11) 
J 

where  Pso  and  Dso  are,  respectively,  the  peak  overpressure  and 
the  duration  of  the  free-water  shock  at  the  charge  surface,  and  PgQ 
and  Dg0  are  the  corresponding  values  for  the  free  airhlast.  In 
practice,  further  simplifications  can  he  introduced  because  of  the 
fact  that  ?so  >>  P'  and  Pg0  >>  e^Pg0  >  and  one  ’.rill  get  finally 
from  Equations -9,  10 ,  and  11 
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?  -*/!> 

F  (t)  =  P  a2e  so 
o  so  o- 

-t/D» 

fi  W  (t)  =  m'P'  a  e  s0 

O  O  SO  0 


Substituting  Equation  12  into  Equations  1  and  8,  and  utilizing  the 
free  air  and  free-water  shock  values  at  the  charge  surface  (refer¬ 
ences  11  and  12) ,  one  obtains  the  following  formula  for  the  water 
shock : 


p(X  ,A  ,x)  =  8.1}  x  io3  -f  e"3-9(T"T,)  h(t-t') 
r  z  xd 


(\  \°‘1 

+  1.1  x  io' k  i  e"26r  H(t)  (psi) 


where 


A  =  qW  3  ,  A  =  rW  3  ,  A  =  zW  3  \ft-lb 

Y  Z 


T  =  (t  -  t)W  ^  ,  xf  -  (q/cQ  -  t)W  ^  \ msec- lb 


Regarding  the  parameters  8  and  S  in  the  airblast  formulas  (Equa¬ 
tions  5  and  6) ,  che  value  of  8  can  be  determined  from  the  energy- 
consideration  ($  =  0.87  in  this  case).  The  parameter  S  is  more 
or  less  related  to  the  overall  geometry  of  the  shock  wave,  and  there 
is  no  clear  way  to  estimate  its  magnitude  beyond  the  fact  that  it  is 
small.  Th’is  S  =  0.15  was  found  to  give  a  better  fit  to  the  test 
data,  although  its  variation  is  not  sensitive  to  the  resulting 
pressure  values. 

Those  pressure  formulas  derived  in  the  case  represented  by 
Equations  5,  6,  and  7  with  8  =  0.87  ,  S  =  0.15  for  the  airblast, 
and  Equation  13  for  the  water  shock  are  plotted  in  figs .  3-5  and 
compared  with  various  test  data  (references  1  thmigh  7»  13,  and  lh). 
The  good  agreement  between  the  theoretical  values  and  the  test  data 
in  these  figures  establishes  the  usefulness  of  these  formulas  for 
spherical  HE  charges.  This  agreement  also  indicates  the  overall 
validity  of  the  various  assumptions  introduced  in  the  process  of 
analysis  in  the  previous  section.  This  fact  is  important  since  this 
gives  some  assurance  of  reliability  for  the  general  formulas  given 
by  Equations  5  through  9,  ".hereby  establishing  the  validity  of  the 
formulas  in  the  next  section  for  the  nuclear  explosion  cases ,  for 
which  any  direct  verification  is  not  feasible. 


Point  source  model 


With  particular  attention  to  the  case  of  stronger 
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explosions,  consider  here  their  initial  stages,  as  described  by  a 
similarity  solution,  which  can  be  derived  as  a  result  of  the  follow¬ 
ing  three  major  assumptions. 

First ,  a  point  source  input  is  assumed  in  the  sense  that 
a  finite  amount  of  mechanical  energy  is  released  instantaneously  at 
a  point  on  the  air-water  interface. 

Second,  the  pressures  at  the  shock  fronts  (in  air  and 
water)  are  assumed  to  be  much  higher  than  the  ambient  pressure  in 
the  undisturbed  region ,  so  that  the  latter  quantity  is  negligible , 

Last ,  it  is  necessary  for  the  similarity  relation  to  assume 
that  both  air  and  water  behave  as  polytropic  substances ,  or  in  other 
words,  the  product  of  (pressure)  x  (density)"^  is  a  function  of 
entropy  only,  where  the  exponent  y  takes  different  values  for  dif¬ 
ferent  media.  The  validity  of  using  this  assumption  in  the  airblast 
region  has  been  well  established.  However,  its  use  for  the  water 
environment  in  connection  with  the  second  assumption  above  should  be 
regarded  with  caution.  This  fact  is  important  to  recognize  when 
taking  into  account  the  difference  between  the  initial  stages  in  weak 
and  intense  explosions. 

In  any  case,  the  assumptions  above  lead  to  a  similarity  of 
flow  field  and  reduce  the  basic  system  of  equations  to  a  simpler 
system  by  use  of  similarity  variables ,  and  furthermore ,  the  similar¬ 
ity  makes  it  possible  to  determine  a  deflate  energy  partitioning 
ratio  of  the  airblast  and  water  shock  in  terms  of  the  released 
energy,  which  is  a  constant  in  the  solution  considered  herein. 

By  utilizing  the  small  value  of  e  „  a  successive  approxi¬ 
mation  scheme  for  the  solution  of  the  similarity  equations  was 
developed  and  tlte  system  of  equations  for  the  first  approximation 
solution  was  found  numerically.  This  provides  an  almost  spherically 
symmetric  solution  (reference  15)  for  the  airblast  region  with  its 
distortion  in  the  order  of  e  ,  and  a  very  oblate  water  shock  front, 
which  is  almost  a  plane  wave  parallel  and  close  to  the  intersection 
surface.  This  solution  was  used  to  compute  the  energy  ratio,  which 
gives  8  =  0.79  and  S  =  h  x  10"3  . 

The  input  functions  F0  and  WQ  of  Equation  9  are  found 
(for  this  approximation)  from 

F  (t)  =  P'(R')R'2 
0  s  s  s 


=  Ip;(r;)es 


Since  the  function  Rg(t)  as  well  as  the  pressure  P^(R|)  are 
available  from  the  free  airblast  data,  these  can  be  utilized  here. 

Conclusions 

General  pressure  formulas  for  the  region  not  so  close  to 
the  explosion  source  were  derived  and  their  output  calculations 
were  compared  with  HE  test  data;  the  comparison  showed  very  good 
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agreement  except  for  a  small  region  on  each  side  of  the  interface, 
wherein  both  measurements  and  formulas  appear  less  reliable.  A  point 
source  model  was  also  developed  to  furnish  the  necessary  initial 
conditions  to  the  general  formulas  for  nuclear  explosion  cases. 
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-  LATER  STAGE 

Fig.  1,  Characteristics  of  initial  and  later  stages  resulting 
from  explosion  at  water  surface  (charge  half  submerged) 
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Fig.  2.  Illustration  of  coordinate  system  and  notation 
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Fig.  U.  Reduced  positi 
with  test  data; 
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e  duration  from  Equation  6  compared 
A  =  0,  and  ll»  ft-lh'1//3 
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Fig.  6.  Water  shock  pressure  time  history  from 
Equation  13  compared  with  test  data 


